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Abstract 

In this paper we study the computability of the stable and unstable 
manifolds of a hyperbolic equilibrium point. These manifolds are the 
essential feature which characterizes a hyperbolic system. We show that 
(i) locally these manifolds can be computed, but (ii) globally they cannot 
(though we prove they are semi-computable). We also show that Smale's 
horseshoe, the first example of a hyperbolic invariant set which is neither 
an equilibrium point nor a periodic orbit, is computable. 



1 Introduction 

Dynamical systems are powerful objects which can be found in numerous ap- 
plications. However, this versatility does not come without cost: dynamical 
systems are objects which are inherently very hard to study. 

Recently, digital computers have been successfully used to analyze dynamical 
systems, through the use of numerical simulations. However, with the known 
existence of phenomena like the "butterfly effect" - a small perturbation on 
initial conditions can be exponentially amplified along time - and the fact that 
numerical simulations always involve some truncation errors, the reliability of 
such simulations for providing information about the long-term evolution of the 
systems is questionable. 

For instance, numerical simulations suggested the existence of a "strange" 
attractor for the Lorenz system |Lor63j , and it was widely believed that such an 
attractor existed. However, the formal proof of its existence remained elusive, 
being at the heart of the 14th from the list of 18 problems that the Fields medal- 
ist S. Smale presented for the new millennium [Sma98j . Finally, after a 35 years 



hiatus, a computer-aided formal proof was achieved in |Tuc98j . |Tuc99j . but this 
example painstakingly illustrates the difference between numerical evidence and 
a full formal proof, even if computer-based. 

For the above reasons, it is important to understand which properties can 
be accurately computed with a computer, and those which cannot. We have 
previously dwelled on this subject on our paper j GZa r which focused on "sta- 
ble" dynamical systems having hyperbolic attractors. The latter systems were 
extensively studied in the XXth century and were thought to correspond to the 
class of "meaningful" dynamical systems. This happened for several reasons: 
the "stability" - structural stability - which implies robustness of behavior to 
small perturbations was believed to be of uttermost importance, since it was 
thought that only such systems could exist in nature; there were also results 
(Peixoto's Theorem |Pei62| ) which showed that, in the plane, these systems are 
dense and their invariant sets (which include all attractors) can be fully char- 
acterized (they can only be points or periodic orbits). There was hope that such 
results would generalize for spaces of higher dimensions, but it was shattered by 
S. Smale, who showed that there exist "chaotic" hyperbolic invariant sets which 
are neither a point nor a periodic orbit (Smale's horseshoe [Sma67J) and that 
for dimensions > 3 structurally stable systems are not dense |Sma66j . 

Moreover, the notion of structural stability was shown to be too strong a 
requirement. In particular, the Lorenz attractor, which is embedded in a system 
modeling weather evolution, is not structurally stable, although it is stable to 
perturbations in the parameters defining the system [ViaOOj (in other words, 
robustness may not be needed for all mathematical properties of the system) . 

Despite the fact that systems with hyperbolic attractors are no longer re- 
garded as "the" meaningful class of dynamical systems for spaces of dimension 
n > 3, they still play a central role in the study of dynamical systems. In 
essence, what characterizes a hyperbolic set is the existence at each point of 
an invariant splitting of the tangent space into stable and unstable directions, 
which generate the local stable and unstable manifolds (see Section |2|. 

In this paper we will study the computability of the stable and unstable 
manifolds for hyperbolic equilibrium points. The stable and unstable manifolds 
are constructs which derive from the stable manifold theorem. This theorem 
states that for each hyperbolic equilibrium point xo (see Section [5] for a defi- 
nition), there exists a manifold S such that a trajectory on S will converge to 
xo at an exponential rate as t — > +oo. Moreover every trajectory converging 
to x lies entirely in S. S is call a stable manifold. Similarly one can obtain 
the unstable manifold U, using similar conditions when t —> — oo. (The stable 
and unstable manifolds of a hyperbolic equilibrium point are depicted in Fig. 
[I]) It is important to note that classical proofs on existence of S and U are 
non- constructive. Thus computability (or, for the matter, semi-computability) 
of the stable/unstable manifold does not follow straightforwardly from these 
proofs. 

We address the following basic question: given a dynamical system and 
some hyperbolic equilibrium point, can we compute its stable and unstable 
manifolds? We will show that the answer is twofold: (i) we can compute the 
stable and unstable manifolds S and U given by the stable manifold theorem 
which are, in some sense, local, since there are (in general) trajectories starting 
in points which do not lie in S that will converge to xq (the set of all these point 
united with S would yield the global stable manifold - see Section [2| , but (ii) 
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Figure 1: Dynamics near a hyperbolic equilibrium point. 



the global stable and unstable manifolds are not, in general, computable from 
the description of the system and the hyperbolic point. Since classical proofs 
on existence of S and U are non-constructive, a different approach is needed if 
one wishes to construct an algorithm that computes S and U. Our approach 
aiming at a constructive proof makes use of function-theoretical treatment of 
the resolvents (see the first paragraph of Section 4 for more details). 

A most likely interpretation for these results is that, since the definition of 
hyperbolicity is local, computability also applies locally. However global homo- 
clinic tangles can occur [GH83] . leading to chaotic behavior for such systems. 
So it should not be expected that the global behavior of stable and unstable 
manifolds be globally computable in general, as indeed our results show. 

In the end of the paper we also show that the prototypical example of an 
hyperbolic invariant set which is neither a fixed point nor a periodic orbit - the 
Smale horseshoe - is computable. 

The computability of simple attractors - hyperbolic periodic orbits and equi- 
librium points - was studied in our previous papers |Zho09j . [GZar] for planar 
dynamics, as well as the computability of their respective domains of attraction. 

A number of papers study dynamical systems, while not exactly in the con- 
text used here. For example the papers |Moo90j . [BBKT01) . |Col05j . |BY06| . 
|Hoy07| provide interesting results about the long-term behavior of dynamical 
systems using, among others, symbolic dynamics or a statistical approach. See 
BGZar for a more detailed review of the literature. 



2 Dynamical systems 

In short, a dynamical system is a pair consisting of a state space where the 
action occurs and a function / which defines the evolution of the system along 
time. See HS74 for a precise definition. In general one can consider two kinds 
of dynamical systems: discrete ones, where time is discrete and one obtains the 
evolution of the system by iterating the map /, and continuous ones, where the 
evolution of the system along time is governed by a differential equation of the 
type 

x = /(*) (1) 
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where t is the independent variable (the "time") and x denotes the derivative 
dx(t)/dt. Continuous-time systems can be translated to discrete-time using 
time-one maps and, in some cases, the Poincare map, and vice-versa (using the 
suspension method). Therefore to study dynamical systems one can focus on 
continuous-time ones. 

Although the essential feature of hyperbolic attractors is the existence of 
an invariant splitting of the tangent space into stable and unstable directions 
generating the local stable and unstable manifolds, for a hyperbolic equilibrium 
point Xo, it can be described equivalently in terms of the linearization of the 
flow around Xq- We recall that Xq is an equilibrium point of ([T]) iff J(xq) = 0. 

Definition 1 An equilibrium point Xq of is hyperbolic if none of the eigen- 
values of Df(xo) has zero real part. 

According to the value of the real part of these eigenvalues, one can determine 
the behavior of the linearized flow near xq: if the eigenvalue has negative real 
part, then the flow will converge to xq when it follows the direction given by 
the eigenvector associated to this eigenvalue. A similar behavior will happen 
when the eigenvalue has positive real part, with the difference that convergence 
happens when t —> — oo. See Fig. [T] The space generated by the eigenvectors 
associated to eigenvalues of Df(xo) with negative real part is called the stable 
subspace E s D ^, Xg y and the space generated by the eigenvectors associated to 
eigenvalues with positive real part is called the unstable subspace Epf,,. 

We now state the Stable Manifold Theorem (as it appears in |Per01j ). 

Proposition 2 (Stable Manifold Theorem) Let E be an open subset of R™ 
containing the origin, let f G C 1 (£'), and let cf> t be the flow of the system |Ip. 
Suppose that /(0) = and that Df(0) has k eigenvalues with negative real part 
and n — k eigenvalues with positive real part (i.e. is a hyperbolic equilibrium 
point). Then there exists a k-dimensional differ entiable manifold S tangent to 
the stable subspace E S D ^^ such that for all t > 0, (f>t(S) C S and for all xq G S 

lim (f> t (x ) = 0; 

t— ^+oo 

and there exists an n ~ k dimensional differ entiable manifold U tangent to the 
unstable subspace E^t, Q s such that for all t < 0, 4>t{S) C S and for all xq G U 



lim 4> t {x ) = 0. 



From the local stable and unstable manifolds given by the Stable Manifold 
Theorem, one can define the global stable and unstable manifolds of ([T]) at a 
hyperbolic equilibrium point x$ by 

oo 

W'jixo) = |J cf-AS) (2) 

3=0 

OO 

Wf(x ) = Q ^(U), 

3=0 

respectively, where 

4> t (A) = {x(t)\x is a solution of with x(0) G A}. 
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We note that W s {x ) and W u (x ) are F a -sets of K™ (recall that a subset F C R" 
is called an i^-set if F = Ujlo ^i' wnere Aj, j £ N, is a closed subset of K n ). 

We end this section by introducing Smale's horseshoe |Sma67| . We refer the 
reader to |GH83j for a more thorough discussion of this set. In essence, Smale's 
horseshoe appears when we consider a map / defined over S = [0, l] 2 . This map 
is bijective and performs a linear vertical expansion of S, and a linear horizontal 
contraction of S, by factors p > 1 and < A < 1, respectively, followed by a 
folding. 

Definition 3 In the conditions defined above (see \GH83^ or ^HSDU^j for pre- 
cise definitions) , the Smale horseshoe is the set A given by 

a= n 

j = -co 

This set is invariant for the function f (i.e. /(A) = A). 

3 Computable analysis 

Now that we have introduced the main concepts of dynamical systems theory 
we will use, we need to introduce concepts related to computability. The theory 
of computation can be rooted in the seminal work of Turing, Church, and oth- 
ers, which provided a framework in which to achieve computation over discrete 
identities or, equivalently, over the integers. 

However, this definition was not enough to cover computability over continu- 
ous structures, and was then developed by other authors such as Turing himself 
|Tur36) . Grzegorczyk [Grz57j, or Lacombe |Lac55| to originate computable anal- 
ysis. 

The idea underlying computable analysis to compute over a set A is to encode 
each element a of A by a countable sequence of symbols from a finite alphabet 
(called a p-name for a). Each sequence can encode at most one element of A. 
The more elements we have from a sequence encoding a, the more precisely we 
can pinpoint a. From this point of view, it suffices to work only with names 
when performing a computation over A. To compute with names, we use Type-2 
machines, which are similar to Turing machines, but (i) have a read-only tape, 
where the input (i.e. the sequence encoding it) is written; (ii) have a write-only 
output tape, where the head cannot move back and the sequence encoding the 
output is written. For more details the reader is referred to |PER89j . |Ko91j . 
[WeiOOj . 

At any finite amount of time we can halt the computation, and we will have 
a partial result on the output tape. The more time we wait, the more accurate 
this result will be. We now introduce notions of computability over M. n . 

Definition 4 1. A sequence {rk} of rational numbers is called a p-name of 
a real number x if there are three functions a, b and c from N to N such 
that for all k e N 7 r k = (-l) a(fe) jjyifi and 

\r k -x\<^-. (3) 
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2. A double sequence {ri t k}i,keN °f rational numbers is called a p-name for 
a sequence {xi}i e jq of real numbers if there are three functions a, b, c from 
N 2 to N such that, for all k,l e N, r l>k = {-l) a{Lk) c ^ k f +1 and 

1 

\ri, k -xi\ < ^fc. 

3. A real number x (a sequence {x{\i e ^ of real numbers) is called computable 
if it has a computable p-name, i.e. the functions a, b, and c are computable 
or, equivalently, there is a Type-2 machine that computes a p-name of x 
({ x i}i&h respectively) without any input. 

In general, p-names for real numbers do not need to be exactly those de- 
scribed above: different definitions can provide the same set of computable 
points (e.g. the p-name could be a sequence of intervals containing xq with di- 
ameter 1/fc for k > 1, etc.). A notable exception is the decimal expansion of x , 
which cannot be used since it can lead to undesirable behavior |Tur37| because 
this representation does not respect the topology of the space M. 

The notion of the p-name can be extended to points in M" as follows: a 
sequence {{r\ kl r2k, ■ ■ ■ ,?"rifc)}fceN of rational vectors is called a p-name of x = 
(x\,X2, ■ ■ ■ ,x n ) G M" if {rjk}ket$ is a p-name of Xj, 1 < j < n. Using p-names, 
we can define computable functions. 

Definition 5 Let X and Y be two sets, where p-names can be defined for ele- 
ments of X and Y. A function f : X — > Y is computable if there is a Type-2 
machine such that on any p-name of x G X , the machine computes as output a 
p-name of fix) G Y. 

Next we present a notion of computability for open and closed subsets of 
W L (cf. [WciOO , Definition 5.1.15). We implicitly use p-names. For instance, 
to obtain names of open subsets of W 1 , we note that the set of rational balls 
B(a,r) — {x G M" : \x — a\ < r}, where a G Q n and r G Q, is a basis for 
the standard topology over M. n . Thus a sequence {(ak,r k )}keN such that E = 
U^ S(afc, r k ) gives a p-name for the open set E. 

Definition 6 1. An open set E C l n is called recursively enumerable (r.e. 
for short) open if there are computable sequences {a k } and {r k }, a k G Q n 
and r k G Q, such that 

E = U k ° =0 B(a k ,r k ). 

Without loss of generality one can also assume that for any k G N, the 
closure of B(a kl r k ), denoted as B(a k ,r k ), is contained in E. 

2. A closed subset A C R" is called r.e. closed if there exists a computable 
sequence {b k }, b k G Q™, such that {b k } is dense in A. A is called co- 
r.e. closed if its complement A c is r.e. open. A is called computable (or 
recursive) if it is both r.e. and co-r.e. 

3. An open set E C W is called computable (or recursive) if E is r.e. open 
and its complement E c is r.e. closed. 

4-. A compact set K C M" is called computable if it is computable as a closed 
set and, in addition, there is a rational number b such that \\x\\ < b for 
all x G K. 
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In the rest of the paper, we will work exclusively with C 1 functions / : E — > 
W 1 , where E is an open subset of M. n . Thus it is desirable to present an explicit 
p-name for such functions. 

Definition 7 Let E = 1J^L B(a,k, 77.), <z& € Q n and r/. £ Q, 6e an open subset 
of W 1 (assuming that the closure of B{au,rk) is contained in E) and let f : 
E — > R" be a continuously differ entiable function. Then a (C 1 ) p-name of f is 
a sequence {Pi} of polynomials (Pi : M. n —¥ R"J with rational coefficients such 
that 

d C i(E)(f,Pi) <2~ l for allleN 

where 

dc HE) (f,Pi)-2_,2 { 1 + llf „ mk + 1 + llDf _ Dmk 



\\g\\k= max \g(x)\. 

x£B(a k ,r k ) 

We observe that this p-name of / contains information on both / and Df in 
the sense that (Pi, P2, ■ . .) is a p-name of / while (DP\, DP2, . . .) is a p-name 
of Df. See |ZW03j for further details. 

Throughout the rest of this paper, unless otherwise mentioned, we will as- 
sume that, in ([lj, / is continuously differentiable on an open subset of W and 
we will use the above p-name for /. 



4 Computable stable manifold theorem 

The stable manifold theorem states that near a hyperbolic equilibrium point xq, 
the nonlinear system 

x = f(x(t)) (4) 

has stable and unstable manifolds S and U tangent to the stable and unstable 
subspaccs and of the linear system 

x = Ax (5) 

where A = Df(x ) is the gradient matrix of / at Xo- The classical proof of the 
theorem relies on the Jordan canonical form of A. To reduce A to its Jordan 
form, one needs to find a basis of generalized eigenvectors. Since the process of 
finding eigenvectors from corresponding eigenvalues is not continuous in general, 
it is a non-computable process. Thus if one wishes to construct an algorithm 
that computes some S and U of @ at Xq, a different method is needed. We 
will make use of an analytic, rather than algebraic, approach to the eigenvalue 
problem that allows us to compute S and U without calling for eigenvectors. 
The analytic approach is based on function-theoretical treatment of the resol- 
vents (see, e.g., [SN42] . |Kat49j . |Kat50j . and |Rob95j 1. 



Let us first show that the stable and unstable subspaces are computable from 
A for the linear hyperbolic systems x = Ax. We begin with several definitions. 
Let 2l# denote the set of all n x n matrices such that the linear differential 
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equation x — Ax, x £ R n , is hyperbolic, where a linear differential equation x — 
Ax is defined to be hyperbolic if all the eigenvalues of A have nonzero real part. 
The Hilbert-Schmidt norm is used for A £ 2l H : ||A|| = Jf^=i E"=i M 2 > 

where a,j is the ij^ n entry of A. For each A £ 2t#, define the stable subspace 
E^ and unstable subspace E^ to be 

E^ = span {v £ R" : v is a generalized eigenvector for an eigenvalue A of A 
with Re(X) < 0} 

E^ = span {v £ R™ : v is a generalized eigenvector for an eigenvalue A of A 
with Re(X) > 0} 

Then R n = E^0E^. The stable subspace E S A is the set of all vectors which 
contract exponentially forward in time while the unstable subspace E^ is the 
set of all vectors which contract backward in time. 

As mentioned above, the process of finding eigenvectors from corresponding 
eigenvalues is not computable; thus the algebraic approach to E^ and E^ is 
a non-computable process. Of course, this doesn't necessarily imply that it is 
impossible to compute E^ and E^ from A, but rather this particular classical 
approach fails to be computable. So even for the linear hyperbolic system x = 
Ax, one needs a different approach to treat the stable/ unstable subspace when 
computability concerned. The approach used below to treat the stable/unstable 
subspace is analytic, rather than algebraic. Let pa (A) be the characteristic 
polynomial for A, 71 be any closed curve in the left half of the complex plane 
that surrounds (in its interior) all eigenvalues of A with negative real part and 
is oriented counterclockwise, and 72 any closed curve in the right half of the 
complex plane that surrounds all eigenvalues of A with positive real part, again 
with counterclockwise orientation. Then 

P A R n = E S A , P%R n = E A 

where 

P\v = -L / (£/ - A)~ 1 vd£, P\v = -L / (£/ - A)- l vd£, 
(See Section 4.6 of [Rob95].) 

Theorem 8 The map H s : 2l# -» A(R n ) = {X\X C R" is a closed subset ofW 1 }, 
A 1 y E S A , is computable, where 2ly is represented " entry wis e" : p = {pij^j—i is 
a name of A if pij is a p-name of aij, the ijth entry of A. 

Proof. First we observe that the map A i-> the eigenvalues of A, A £ 2l#, 
is computable. Assume that Ai, . . . , Xk, Pk+i, ■ ■ ■ , Pn (counting multiplicity) are 
eigenvalues of A with Re(Xj) < for 1 < j < k and Re(pj) > for k+1 < j < n. 
Then from the p-names of the eigenvalues, one can compute two rectangular 
closed curves ^y A and j A , where 7^4 is in the left half of the complex plane that 
surrounds all eigenvalues Xj with 1 < j < k and 7^ is in the right half of the 
complex plane that surrounds all eigenvalues pj with k + 1 < j < n, and both 
j A and 7^ are oriented counterclockwise. From j A and 7^ one can further 
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computes the maps P\,P\ : R n -> K n , where P\v = A- jy (£1 - A)- X vd^ 
ju = ^jj J^ 2 (£7 — ^4) _1 wd^. We note that, on the one hand, E^ = 



(P|) _1 ({0}): and on tne otner nand > E A = e a = P\ Rn i where if denotes the 
closure of the set K. Then by Theorem 6.2.4 of |Wei 00j . W A is both r.e. and co- 
r.e. closed, thus computable. The same argument shows that EJ4 is computable 
from A. m 

Next we present an effective version of the stable manifold theorem. Consider 
the nonlinear system 

x = /(*(*)) (6) 

Assume that ^ defines a dynamical system, that is, the solution x(t,Xo) to (6) 
with the initial condition x(0) = Xq is defined for all ieM. Since the system (6) 
is autonomous, if p is a hyperbolic equilibrium point, without loss of generality, 
we may assume that p is the origin 0. 



Theorem 9 Let f : R n — > R n be a C 1 -computable function (meaning that both 
f and Df are computable). Assume that the origin is a hyperbolic equilibrium 
point of such that Df(0) has k eigenvalues with negative real part and n — k 
eigenvalues with positive real part (counting multiplicity), < k < n. Let 
x(t, x ) denote the solution of |6p with the initial value x at t — 0. Then 
there is a (Turing) algorithm that computes a k-dimensional manifold S C K.™ 
containing such that 

(i) For all Xq £ S, lim t _j. +00 x(t, Xq) = 0; 

(ii) There are three positive rational numbers 7, t, and 8 such that \x{t, Xq)\ < 
r )2~ et for all t > whenever x$ € S and \xq\ < 5. 

Moreover, if k < n, then a rational number rj and a ball D with center at the 
origin can be computed from f such that for any solution x(t,Xo) to the equation 
@ with x eD\S, {x(t, x ): t > 0} (£ B(0, 77) no matter how close the initial 
value xq is to the origin. 

Proof. First we note that under the assumption that / is C 1 -computable, the 
solution map x : R x E™ -> R n , (t, a) ^ x(t, a), is computable ( [GZB09j ). 

Let f{x) =Ax + F(x), where A = Df(0) and F(x) = f{x) - Ax. Then the 
equation (JoJ) can be written in the form of 

x = Ax + F(x) (7) 

The first step in constructing the desired algorithm is to break the flow e At gov- 
erned by the linear equation x = Ax into the stable and unstable components, 
denoted as Ir 1 (t) and Jp 2 (t) respectively. By making use of an integral formula, 
we are able to show that the breaking process is computable from A. The details 
for the first step: Since F(x) = f(x) - Ax, it follows that F(0) = 0, DF(0) = 0, 
F and DF are both computable because / and Df are computable functions 
by assumption. Thus there is a computable modulus of continuity d : N — > N 
such that 

|-F(s:) - F(y)\ < 2~ m \x - y\ whenever |x| < 2- d ^ & |y| < 2~ d ^ (8) 
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Since Df is computable, all entries in the matrix A are computable; conse- 
quently, the coefficients of the characteristic polynomial det(A — XI n ) of A are 
computable, where det(A— XI n ) denotes the determinant of A— XI n and is the 
n x n unit matrix. Thus all eigenvalues of A are computable, for they are zeros of 
the computable polynomial det(^4 — A/„) . Assume that Ai , . . . , A& , /ifc+i, ... ,fx n 
(counting multiplicity) are eigenvalues of A with Re(Xj) < for 1 < j < k and 
Re([ij) > for k + 1 < j < n, where Re(z) denotes the real part of the complex 
number z. Then two rational numbers a > and a > can be computed from 
the eigenvalues of A such that Re(Xj) < — (a + cr) for 1 < j < k and Re(/ij) > a 
for k + 1 < j < n. We break a into two parts for later use: Let ot\ and a 2 be 
two rational numbers such that 

< ai < a and a\ + a 2 = a. (9) 

Let M be a natural number such that M > max{a+cr, 1} and max{|Aj|, \/ii\ : 
I < j < k, k + 1 < I < n} < M — 1. We now construct two simple piecewise- 
smooth close curves Li and T 2 in IR 2 : Ti is the boundary of the rectangle with 
the vertices (-a - a, M), (-M, M), (-M, —M), and (-a - a, -M), while L 2 
is the boundary of the rectangle with the vertices (cr, M), (M,M), (M, — M), 
and (cr, —M). Then Li with positive direction (counterclockwise) encloses in its 
interior all the Xj for 1 < j < k and T 2 with positive direction encloses all the /ij 
for k + 1 < j < n in its interior. We observe that for any £ £ Li IJ T 2 , the matrix 
A-£I n is invertible. Since the function g : Fx \JT 2 ->• E, g(£) = | j (^4 — ^J„) _1 | j, 
is computable (see for example [Zho09 ), where (^4 — £/ n ) _1 is the inverse of the 
matrix A — £/„, the maximum of g on Ti 1JL 2 is computable. Let K\ € N be 
an upper bound of this computable maximum. Now for any feR. from (5.47) 
of |Kat95j . 



~\ e t-XA-tl n )-^-~\ e t\A-tl n )-\ 
2m J Tl 2m J V2 

IrAt) + hS) 



We recall that e is the solution to the linear equation x — Ax. Since A is 
computable and integration is a computable operator, it follows that Ir t and Ir 2 
are computable. A simple calculation shows that || — e t ^(A—^I n )^ 1 d£ l \\ < 

iK^Ie-^+^/n for t > and || - f F2 e^(A - £/„) _1 d£|| < AKxMe^/iT 
for t < 0. Let 

K = 4MK t (11) 

Then 

1 1%(*)| I < Ke- (a+a '> t for t > and ||Ir 2 (*)|| < Ke at for t < (12) 



The two estimates in ( 12 ) show that I Fl (t) and Ir 2 (t) are stable component and 
unstable component of e A . The first step is now complete. 

The second step in the construction is to compute a ball, B(0, r), surrounding 
the hyperbolic equilibrium point such that 5(0, r) contains a set of potential 
solutions x(t, xo) to ([7| satisfying the conditions (i) and (ii) described in Theo- 
rem [9} The details for step 2: 
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Claim 1. Consider the integral equation 



u(t, a) = Jr/i (t)a + / Ir 1 (t — s)F(u(s,a))ds ■ 



Ir 2 (t-s)F(u(s,a))ds (13) 



where the constant vector a is a parameter. If u(t, a) is a continuous solution 
to the integral equation, then it satisfies the differential equation Q with the 
initial condition u(0, a) — Ir 1 (0)a — Ir 2 (—s)F(u(s, a))ds. 

Proof of Claim 1. See Appendix 1. 

We observe that the solution to the integral equation is the fixed point of the 
operator defined on the right hand side of the equation (13 1. Next we compute 
an integer mo such that 2~ m ° < j^, and then set 



r = 2- d{rno) /2K, B = B(0, r) = {a G 
where the computable function d is as in (Jsl). 



\a\ < r} 



(14) 



Claim 2. For any a £ B and t > 0, define the successive approximations as 
follows: 



u(°\t,a) =0 

u {j) (t, a) = Z ri [t)a + [ Z ri (t - s)F{u ij -V ( s , a ))ds 
Jo 

Ir 2 (t- S )F(u^ 1 \s,a))ds, j>l 



(15) 



then the following three inequalities hold for all j G N: 

\u U) (t, a) - u {j - 1] (t, a)\ < K\a\e~ ait 
\u {i \t,a)\ < 2- d ^e' ait 

and for any a G B, 

\u^(t,a)-u (j) (t,a)\ < 3A"|o-o| 
Proof of Claim 2. See Appendix 2. 



(16) 
(17) 

(18) 



It then follows from ( 15 1 and ( 16 ) that {u^ (t, a)}°^ 1 is a computable Cauchy 
sequence effectively convergent to the solution u(t, a) of the integral equation 
( 13 1, uniformly in t > and a G B. Consequently the solution, t, a i-> u(t, a), is 
computable. Furthermore, (17) and (18) imply that for alH > and a, a G B, 



\u(t, a)\ < 2- di - mo) e - ai \ \u(t, a) - «(t, S| < 3K\a - a\ 



(19) 



Thus limt-^oo u(t, a) = for all a G B. Moreover, the first inequality in (19) 
shows that u(t, a) satisfies condition (ii) of Theorem [9] with 
and <5 = 2- d ^/2K. 



Although Claim 2 shows that for any a G B, the integral equation (13) 
has a computable solution u(t, a) and this solution satisfies the conditions (i) 
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and (ii) of Theorem [9j we may not take B as a desired stable manifold S be- 
cause if u(0, a) 7^ a, then u(t,a) is not a solution to the equation Q with the 
initial value a at £ = 0. Nevertheless, the set P> provides a pool of potential 
solutions to ^ on a stable manifold. This leads us to the next step of the proof. 

The last step in the construction of the desired algorithm is to extract a 
fc-dimensional manifold S from B, using a computable process, such that for 
any a g S, u(0, a) = a. Then if we set x(t, a) = u{t, a) for a G 5, by Claims 
1 and 2, x(t, a) is the solution to the differential equation Q with the initial 
value a at t = and satisfies the conditions (i) and (ii) of Theorem [9] Now for 
the details. 

Let P 1 = J ri (0) = -gij / ri (A - ^r 1 ^ and P 2 = I Ta (0) = -£tf Ta (A - 
£,I n )- l d£. Then P x R n = E^, P 2 M" = E%, 1" =E^0E^ with dimE^ = k and 
dimE^ — n — k, PjPk = SjkPj {5jk = 1 i£ j = k and Sjk = if j ^ k), and 
P 1 + P 2 = I is the identity map on R n (c.f. §1.5.3 and §1.5.4, [Kat95j : §4.6 of 
|Rob95| ). Since A is computable, so are Pi and P 2 . Moreover, I ri (t)P 2 = 
for any < £ t as the following calculation shows: Let R((,) denote (A — ^/„) _1 . 
Then we have 

- = - e 2 /«)H(6) - r^i){a - fi/„)ij(e a ) 

= fife) (6-6)^(6) 
= (6-6)fi(6)fi(6)- 



Using the last equation we obtain that for any v € R™, 

i ri (t)P2v= (J^j jf e «*P(0^^ 



(20) 



1 ^ 2 



2™ 



ri Jr 



27TI y 



2 

v 27Ti / 

A similar computation shows that for any igi, 

7r 2 (i)Pi« = and P 2 Ir 2 {t)v = h 2 (t)v, v € R n (21) 
Now let us use these results to compute the projections of u(0, a) in E^ and E^: 
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for any a G K. n , 

Piu(0,a) =Pi (l Tl (0)a- J Ir 2 (s)F(u(s,a))ds 



and 



/>oo 

PiPia- / P 1 Ir 2 (~s)F(u(s,a))ds 
Jo 

Pi a (22) 



P 2 u(0,a) =P 2 f/ ri (0)o- J Ir 2 (-s)F(u(s,a))ds 
= P 2 P ia -P 2 n J r3 (-s)P(u(s,a))ds 



Jr 2 (~ s )P( m ( s j a))ds 



We note that for any a G 



I Tx {t)a = I Vl {t){Pia + P 2 a) 

= I r ,{t)Pxa + I Tl {t)P2a 
= I Vl {t)P x a 



which implies that if the solution u(t, a) of the integral equation ( 13 1 is con- 
structed by successive approximations (15), then P 2 a does not enter the com- 
putation for u(t, a) and thus may be taken as zero. Therefore the projection of 
w(0, a) in E^ satisfies the equation 

/•oo 

P 2 u(Q,a) = - Ir 2 {-s)F(u(s,P 1 a))ds (23) 
Jo 

Next we define a map 4>a ■ E s A {r) — > E A , b H> — J °° Ir 2 (—s)F(u(s,b))ds for 
b G E s A (r), where r is the rational number defined in ( [l4| and E s A (r) = {b G 

: |6| < r/2}. We observe that the compact set E A (r) is computable since 
the closed set E^ is computable (proved in Theorem [8]). Obviously the map 
4>A is computable. By Theorem 6.2.4 |WeiOO| the compact set (^[E^r)] is 
computable. 

Now we are ready to define S: 

S = {b + Mb)-beE s A (r)} 

Then 5* is a manifold of dimension k. For any a G S, a = b + </>a(&) for some 
b G E s A (r). Since R n = E^©E^, b G E^ and <p A (b) G E A , it follows that 

6 = Pia and A (6) = ^(Piflt) = P 2 a (24) 



Combining (|22j, (J23J), and (24) we obtain that for any a £ S, u(0,a) — 
Piu(0, a) + P 2 u(0, a) = P x a + 4>a{Pio) = Pia + P 2 a = a. Thus for any a G S, 
u(0, a) = a; that is, a is the initial condition of u(t, a) at £ = 0. 

Claim 3. S is a computable closed subset of M. n . 
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Proof. Since E^(r) is computable, there is a computable sequence {bj} C E s A (r) 
that is effectively dense in K s A (r); that is, there is a computable function i/j : N — > 
N such that E^(r) C 2~ fe ) for all jfe e N (c./. |Zho96j V The following 

estimate shows that the computable sequence {bj + 4>A{bj)} is effectively dense 
in S, thus the closed manifold S is computable. For any 6+</>^(6), b + 4>A(b) G S, 
we have 

\b + Mb)-(b + Mb))\ 
< 



< 



< 



Ir 2 (-s)F{u{s,b))ds+ / Ir 2 {s)F(u(s,b))ds 
o Jo 



\\Ir 2 (-s)\\-\F( U (s,b))-F(u( S ,b))\ds 



b\ + \Mb)-Mb)\ 

b\ 

b\ 
b\ 
b\ 
b\ 



o 

2 -m | 

if 
cr2 m o 

is: 

<t2 to ° 



m(s,6) -m(s,6)| / Ke- as ds 
Jo 

\u(s,b) - u(s,b)\ 



3K\b-b\ 



The estimates (U2|), (|8|), (14), and (18) are used in the above calculation. The 



proof of claim 3 is complete. 

For every a <E S, set x(t,a) — u(t,a); then by Claims 1, 2, and 3, x{t,a) 
is the solution to the differential equation ^ with initial value a at t = and 
x{t,a) satisfies the conditions (i) and (ii) of Theorem[9] 

Finally we show that if x(t,xo) is a solution to the differential equation ^ 
satisfying < |x | < 2~ d( - m ^/AK 2 but x g S, then there exists t' > such 
that |x(t',xo)| > 2~ d ( m °K This proves the last part of the theorem if we set 
rj = 2- d ( m °) and D = {x e E" : |jc| < 2" d ( m °) /4K 2 }. Indeed, if otherwise 
1^(^,^0)1 < 2 _d ' m °) for all t > 0. We show in the following that this condition 
implies xo € 5, which is a contradiction. 

Since x(i,Xo) is the solution to x = Ax + F(x) with the initial value Xo, it 
follows that (c./. Theorem 4.8.2 |Rob95j ) 



x(t,x ) 



e At x 



,{t-s 



)A F(x{s,x ))ds 



which, using (10), can be rewritten as 
x(t,x ) 



(25) 



Iri(t)xo +I T2 {t)x + I I ri (t- s)F(x(s,x ))ds+ / Ipjt - s)F(x(s,x ))ds 



Ir 2 (t - s)F(x{s, x„))ds + h 2 (t)b 



= I ri (t)x + / I^it- s)F(x(s,x ))ds 



where b = x + /„ I r2 (-s)F(x(s,x ))ds (c.f. §1.5.3 |Kat95j ). Note that b 
is well defined since we assume that \x(t, x )| < 2~ d ( m °) for all t > 0, then, 



from (|8]), \F(x(t,Xo))\ is bounded for all t > 0; consequently, by (12), the 
integral Jr 2 (-s)F(x(s, xo))ds converges. We also note that the first three 
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terms in the above representation for x(t, xq) are bounded. Moreover, Ir 2 (t)b 



Ir 2 (t)P 1 b + h 2 {t)P 2 b = Ir 2 {t)P 2 b since Ir 2 (t)P 1 b = by ([21]). We claim that 
if P 2 b ^ 0, then Ir 2 (t)b is unbounded as t — > oo. We make use of the residue 
formula to prove the claim. Recall that T 2 is a closed curve in the right-hand 
side of the complex plane with counterclockwise orientation that contains in its 
interior fij, k + 1 < j < n, where fij are the eigenvalues of A with Re([ij) > 0, 
which are the exact singularities of R(£) = (A — ^In) 1 in the right complex 
plane. Then by the residue formula, 

1 f " 

— e*R(0d£ = E eHt ™ s H R (26) 

m ^ V2 l=k+l 

where res^.i? is the residue of R at Since Re(/j,j) > 0, if P 2 b ^ 0, then 
Jr 2 (t)6 = /r 2 (t)P 2 6 



~ ( e t HA-ZI n )- 1 P 2 bdt 



= _ J2 e^ t res H RP 2 b 

l=k+l 

is unbounded as t — > oo. This is however impossible because the first three terms 



in (25) are bounded and we have assumed that \x{t, Xq)\ < 2 d ( m °) for all t > 0. 



Therefore, P2& = 0; consequently, Ir 2 (t)b — Ir 2 {t)P 2 b = 0. This last equation 



together with (25 1 shows that x(t,Xo) satisfies the integral equation (13). Now 



let x'(t, Xq) be the solution to the integral equation (|13[) with parameter xq and 



constructed by the successive approximations (15). Then 



x'(t,x )=I ri {t)x + / Ir 1 (ts)F(x'{s,x ))ds- / Ir 2 {t-s)F{x'(s,P 1 x ))ds 



By the uniqueness of the solution, x(t,Xo) = x'(t,xo); in particular, xq — 
x(0,Xq) = x'(0,xq) = P1X0 — Ir 2 {~s)F(x' (s, PiXo))ds. Since P\Xq € 
and \\PiX \\ < ||i\||||aso|| < K ■ 2 d ^/4K 2 = r/2 (recall that \\Pi\\ < K by 



(12)), it follows that P\x$ g W A (r), which further implies that xq = P\Xq — 
Ir 2 (—s)F(x' (s, P\Xo))ds = P\Xq + ^>a(P\Xq) € S. This contradicts the fact 
that Xq is not on S. The proof is complete. ■ 

Theorem 10 Let f : ffi™ — > M n 6e a C 1 -computable function. Assume that the 
origin is a hyperbolic equilibrium point of |6|) such that Df(0) has k eigenvalues 
with negative real part and n — k eigenvalues with positive real part ( counting 
multiplicity), < k < n. Let x(t,xo) denote the solution to |#|) with the initial 
value xq at t = 0. Then there is a (Turing) algorithm that computes a (n — k)- 
dimensional manifold U C R™ containing such that 

(i) For all x a £ U, limj^-oo x(t,x ) = 0; 

(ii) There are three positive rational numbers 7, e, and S such that \x(t,xo)\ < 
72 e * for all t < whenever xq € U and \xq\ < 5. 

Moreover, if k > 0, then a rational number 77 and a ball D can be computed 
from f such that for any solution x(t,Xo) to the equation with x a G D\U, 
{x(t,xo) : t < 0} (jL 5(0,77) no matter how close the initial value xq is to the 
origin. 
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Proof. The unstable manifold U can be computed by the same procedure as 
the construction of S by considering the equation 

x = -Ax - F(x(t)) 

m 

The proof can be easily extended to show that the map: §h — > J? x .ft, 
/ H> (Sf, Uf), is computable, where $h is the set of all C 1 functions having the 
origin as a hyperbolic equilibrium point, A is the set of all compact subsets of 
R™, and Sf and Uf are some local stable and unstable manifolds of / at the 
origin respectively. 

5 Computability and non-computability of global 
stable/unstable manifolds 

Although the stable and unstable manifolds can be computed locally as shown 
in the previous section, globally they may not be necessarily computable. 

In |Zho09| it is shown that there exists a C°° and polynomial-time com- 
putable function / : E 2 — > E 2 such that the equation x = f(x) has a sink at 
the origin and the basin of attraction (also called the domain of attraction) , 
Bf(0), of / at is a non-computable open subset of E 2 . Since a sink is a hy- 
perbolic equilibrium point (with all eigenvalues of the gradient matrix 13/(0) 
having negative real parts) and the basin of attraction at a sink is exactly the 
global stable manifold at this equilibrium point, we conclude immediately that 
the global stable manifold is not necessarily computable. 

One the other hand, it is also shown in |Zho09j that, although not necessarily 
computable, the basin of attraction Bf(xo) is r.e. open for a computable sink 
x Q of a computable system /. In other words, the open subset Bf(x ) can 
be plotted from inside by a computable sequence {B n } of rational open balls, 
Bf(xo) = (J B n , but one may not know the rate at which Bf(xo) is being filled 
up by -B n 's if Bf(xo) is non-computable. When an equilibrium point xq of the 
system x = f(x) is a saddle (not a sink nor a source), then the global stable 
(unstable) manifold Wi(xo) (Wf(xo)) is in general an F CT -set of E n . As in the 
special case of the basin of attraction, Wi (xq) can also be plotted from inside by 
a computable process. To make the concept precise, we introduce the following 
definition. Let .F(R") denote the set of all F CT -subsets of E" . 

Definition 11 A function f : X — > is called semi-computable if there is 

a Type-2 machine such that on any p-name of x G X, the machine computes as 
output a sequence {aj t k}, 0>j,k G Q" > such that 

oo 

f(x) = |J K* : k E N} 

3=0 

where A denotes the closure of the set A. 

We call this function semi-computable because we can plot a dense subset 
of the set f{x) by a computable process, but we cannot tell in a finite amount 
of time what the "density" is. 
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Theorem 12 The map ip : ^(R^xl" -> T(R n ), (f,x ) -» Wf(x ), is semi- 
computable, where (f,Xo) G dom(?p) if Xq is a hyperbolic equilibrium point of f , 
and Wi(xq) is the global stable manifold of f at xq. 

Proof. From Theorem [9] one can compute from / and Xo a compact subset S 
of R n , which is a local stable manifold of / at xq. Note that the global stable 
manifold of / at Xq is the union of the backward flows of S, i.e., 

oo 

Wf{x ) = (J <t>-i{S) 
j=o 

where 4>t{o) is the flow induced by the equation x = fix) at time t with the 
initial data x(0) = a. Since the sequence {ip-j(a)} is computable from / and a 
GZB09 , it follows that the sequence {(f>-j(S)} of compact sets are computable 
from / and xq (Theorem 6.2.4 |WeiOO| ). In particular, a sequence {aj,k}j,keH C 
R n can be computed such that Wj(xq) — Ujlo { a j-k '■ k € N}. ■ 

The function / in the counterexample mentioned at the beginning of this 
section is C°° but not analytic. It is an open problem whether or not the 
stable/unstable manifold(s) of a computable analytic hyperbolic system x — 
f(x) (i.e. f is computable and analytic) is computable. In the following we 
present a negative answer to a weaker version of the open problem; we show that 
it is impossible to compute uniformly the closure of the global stable/unstable 
manifold from (/, xq), where / is analytic and xq is a hyperbolic equilibrium 
point of /. Let's denote by the set of real analytic functions / : K™ — > E n 

and A(M. n ) the set of all closed subsets of R™. With Wijsman topology W on 
A(R n ), (A(M. n ),W) is a topological space that is separable and metrizable with 
a complete metric. It can be shown that the Wijsman topology on A(M. n ) is 
the same as the topology induced by the p-names giving rise to the notion of 
computable close subsets of R n as defined in Def. [6] ([BW99 ). Thus if a map 
F : X — > A(M. n ) is computable (with respect to above p-names for elements of 
.4.(11")), then F is continuous (with respect to Wijsman topology on A(M. n )) 
(Corollary 3.2.12 of [WeiOOj ). 

Theorem 13 The map tp : u(R n ) x R" -> A(R n ), (f,x ) -> Wf(x Q ) (the clo- 
sure ofWJ(xo)), is not computable, where (f,Xo) € dom(tp) ifxo is a hyperbolic 
equilibrium point of f , and WJ(xq) is the global unstable manifold of f at Xq. 

Proof. Consider the following system x — f^x) taken from |HSD04| . 

x = x 2 - 1 (27) 
?/ = -xy + n{x 2 - 1) 



The system (271 has two equilibria: z\ = (—1,0) and Z2 = (1,0). Df^(zi) has 
eigenvalues —2 and 1, associated to the eigenvectors (—1,-2/^/3) and (0,1), 
respectively, and Df^(zi) has eigenvalues 2 and —1, associated to the eigenvec- 
tors (—1,-2^/3) and (0,1), respectively. Therefore both points Z\ and z 2 are 
saddles (the behavior of the system is sketched in Fig. [2]). From this information 
and the fact that any point (— l,y) or (l,y) can only move along the y axis, 
one concludes that the line x — — 1 is the unstable manifold of Z\ and the line 
x = 1 gives the stable manifold of z<i. Note that z\, z-i and the above manifolds 
do not depend on fi. 



17 



Let us now study the unstable manifold of 22- It is split by the stable 
manifold x = 1, and hence there is a "right" as well as a "left" portion of 
the unstable manifold. Let us focus our attention to the "left" portion. From 



the system (27) one concludes that any point z near Z2, located to its left (i.e. 
z < Z2) will be pushed to the line x = —1 at a rate that is independent of the 
y-coordinate of z. 

When fj, = 0, the eigenvector of Df fl (z2) associated to the unstable manifold 



is (—1,0). Looking at (27), one concludes that the "left" portion of the unstable 



manifold of Z2 can only be the segment of line 

Ulo = {(x,0) e K 2 | - 1< x < 1}. 

Now let us analyze the case where \i < 0. Since (—1,-2^/3) is the eigen- 
vector of .D/ M (z2) associated to the unstable manifold, as the unstable manifold 
of Z2 moves to the left, its y-coordinate starts to grow. The "left" portion of 
the unstable manifold is always above the line y = (if it could be y — 0, the 



dynamics of (27) would push the trajectory upwards), and as soon as its IE- 
coordinate is less than (this will eventually happen), the trajectory is pushed 
upwards with y-coordinate converging to +00. Notice that the closer /1 is to 
0, the lesser the y-component of the trajectory grows, and the closer to z\ the 
unstable manifold will be. 

If Ui „ represents the left portion of the unstable manifold of Z2 and 

A = U l>0 U{(-l,y) eM 2 |0<y} 

then one concludes that 

lim d(A, Uia) = 

iu->-0~ 

where d is the Hausdorff distance on R". 

Suppose that i\> is computable. Then, in particular, the map x '■ K — >• A(M. n ) 
defined by x(/x) = WJ (22), is also computable. Since computable maps must 
be continuous (Corollary 3.2.12 of [WeiOOj ) . this implies that \ should be a 
continuous map. But the point /1 = is a point of discontinuity for \ since 
(without loss of generality, we restrict ourselves to the semi-plane x < 1) 

lim x(m) = lim U Lfi = A^ U L0 = x(0). 

fi — y fj,— >0~ 

and hence this map cannot be computable, which implies that ip is not com- 
putable either. ■ 



6 The Smale horseshoe is computable 

In this section we show that Smale's horseshoe is computable. 

Theorem 14 The Smale Horseshoe A is a computable (recursive) closed subset 
in I = [0,1] x [0,1]. 

Proof. We show that f2 = I\ A is a computable open subset in / by making use 
of the fact |Zho96| : An open subset U C I is computable if and only if there is a 
computable sequence of rational open rectangles (having rational corner points) 
in /, { Jk}kL , such that 
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[i < n = o 

Figure 2: Dynamics of the system (27). 



M > 



(a) J k cU for all fc e N, 

(b) the closure of Jk, Jf., is contained in U for all fc € N, and 

(c) there is a recursive function e : N — > N such that the Hausdorff distance 
d(I \ ^k=l J kJ \U)< 2~ n for all n e N. 

Let / : 7 ->■ M 2 be a map such that A = |XL-oo(/"C0 n 7 ) is the Smale 
horseshoe. Without loss of generality assume that / performs a linear vertical 
expansion by a factor of (i = 4 and a linear horizontal contraction by a factor of 

A = j . For each n£N, let t/ n = J \ |X=-„(/ fc 0O n 7 )- Then 7 \ A = U^o U n- 
Moreover, 

(1) there exists a computable function a : N — > N such that U n is a union of 
a(n) rational open rectangles in I, 

(2) rf(A,/\ULo^)<(V4)". 

Define e : N -> N, e(n) = ELo a ( fc )' n e N - Thcn ' lt follows lrom thc 
lemma that I \ A is indeed a computable open subset of I. m 



7 Conclusions 

We have shown that, locally, one can compute thc stable and unstable manifolds 
of some given hyperbolic equilibrium point, though globally these manifolds 
are, in general, only semi-computable. It would be interesting to know if these 
results are only valid to equilibrium points or can be extended e.g. to hyperbolic 
periodic orbits. In [GZarj we provide an example which shows that the global 
stable/unstable manifold cannot be computed for hyperbolic periodic orbits. 
However, the question whether locally these manifolds can be computed remains 
open. 

Appendix 1. Proof of Claim 1. Assume that u(t,a) is a continuous solution 



to the integral equation (13). We show that it satisfies the differential equation 
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Q. To see this, we first establish a relation between Jr^., j = 1,2, and the 
Jordan canonical form of A: 



A = C 



P 
Q 



c- 1 



where C is an n x n invertible matrix, P is a fc x k matrix with eigenvalues Xj , 
1 < j < fc, and Q is a (n — k) x (n — A;) matrix with eigenvalues /ij, + l < j < n. 
By (10), we have 

e M = I Tl {t)+hS) 



On the other hand, it is straightforward to show that 

3 P * 



z M = C 







C' 1 + c 





e«* 



Since the first k columns of C consists of a basis of the stable subspace PiR n 
and the last (n — k) columns of C is a basis of the unstable subspace P 2 M. n , it 
then follows from (201 and (21) that 



e pt 




Thus 



and 



Ir 2 (t)C 
I Fl (t)C 

J ri (t) = C 





e Qt 



C- 1 = 0, and 

c- 1 = 



e pt 




c- 1 



c^ 1 



(recall that R n = PiK" P 2 R n ). Consequently the integral equation ([13| can 
be written in the following form: 



u(t, a) = C 



e pt 




C~ l a- 



C 



e P(ts) o 





C^Fiuis^^ds (28) 



C 





e^K*-*) 



C -1 F(u(s, a))ds 



It is known that if u(t, a) is a continuous solution to (28), then it is the solution 
to (c./. §2.7 [PerOlj ). The proof is complete. 

Appendix 2. Proof of Claim 2. For a, a £ B and t > 0, where £? = -B(0, r) = 
{x e E™ : |x| < r}, r — 2~ d( ~ m °"> /2K (d is the computable function defined in 
and mo is a positive integer such that 2 _m ° < define 

u (0) (i,a) =0 



,0")/ 



(t,a) = Ir 1 { t ) a + / /r 1 (t-s)^(w (j '" 1) (s,a))rfs 



Ir 2 (t- s)F(u (j - 1) (s 1 a))ds, j>l 
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We show that the following three inequalities hold for all j € N: 

\u^(Tt,a)-u^-^(t,a)\ < K\a\e~ ait 

\u^(t,a)\ < 2- d ( m °h~ ait 
|u w (t,a) -u^(t,a)\ < 3K\a-a\ 

We argue by induction on j. Since u^°\t,a) ~ for any a, by (|9| and (12) 
we get \u^{t,a)\ = \\I Tl (t)a\\ < Ke^ 01 ^ 1 ■ 2- d( - ma ~> /2K < 2- d ( m oT e - Qlt , and 
|u (1) (<,a)-" (1) M)l = \I Tl (t)(a-a)\ < Ke- {a+a)t \a-a\ < K\a-a\, the three 
inequalities hold for j = 1. The estimate (12) is used in calculations here. 
Assume that the three inequalities hold for all k < j. Then for k = j + 1, 



u U+1) (t,a) -u (j) {t 


a)\ 


< f\\hAt-s)\\- 

Jo 


\F(u u) (s,a)) - F(u^-^(s,a))|ds 


poo 

+ 1 l|Jr a (t-*)|| 


• \F(u^(s,a)) - Fiu^-^is^Wds 


< f \\IrAt-s)\\- 
Jo 


2- m °\u {j) (s,a)-u < - j - 1 ' ) (s,a)\ds 


poo 

+ J t ||/r a (t-*)|| 


■2- m °\u u) (s,a) -u^-^(s,a)\ds 


< 2-" l ° [ Ke- (ai - 
Jo 


- a2 +a)(t~s) K \a\e~ aiS O _ mo f°° 
2*-i dS + 2 J t 



Ke^ K ^T S ds 

= 2-"i -^ 2 | Q l c -(a 1 +a 2 + CT )t f (aa+o-jg. , g-mo gf f°° {-(ai+<r)sj 

_ 2 - mn ^ 2 kl c - (ai+a , +CT ) t e^ + ^-i ^^ j^r^ 

2J- 1 a 2 + cr 2J- 1 -(ai+cr) 

m if 2 |a| e" Ql * m K 2 |a| e- Qlt 
< 2~ 1 1 [_ 9— mo 1 1 



2J'- 1 a 2 + cr 2J'- 1 ax + a 

< ~^r e ~ ait ( reca11 that 2_mo ^ 5?) 

and furthermore, 
\u {j+1 \t,a)\ 

< \vP(t,a)\ + \u (j+1) (t,a) -u^{t,a)\ 
<|^)(i,a)| + S e -^ 

3 I 

< — — e _aii (induction hypothesis on u 

for k < j) 

fe=i 

< 2K\a\e~ ait < 2K e - ait ■ 2- d{mo) /2K = 2- dt - m ^ e - ait 

Lastly we show that if \u^(t,a) — ?/ fe )(i,a)| < 3K\a — d\ holds for all k < j, 
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then it holds for j + 1. 

\u^ j+1 \t,a) -u^ +1 \t,d) 



- 1 (t)(o-5) + J hAt-s) (F{u^\s,a))-F{u^\s,d)))ds- 
- J Jr 2 (t-«) (F(uW(s,a))-F(MW(s,a)))ds 

< |/ ri (t)(a-a)|+ / ||/ ri (*-«) 
Jo 

/oo 
||Jr 2 (t-«)||- F(«W)(*,a))-F(u^(a,a)) 



F(u U \s,a))-F(u U \s,d)) 
ds 



ds+ 



< K\a-~a\+2- ma 

< K\a- d\+2- m ° -3K\a-a 



u^(a,a)-«^(a,o)| 



^ e -(a+a)(t- s ) ds + jf ^(*-») ( ,. s 



AT 



K 



a + a <t 



= K\a — a\(l + 2~ m ° ™- + 2- m °—) 
\ a + a a ) 

< 3K\a - a\ (Recall that 2- m ° < 

The proof is complete. 
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